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1  13  ABSTRAC  T  1 

The  time-limited  nonnegative  data  windows  that  minimize  the  bias  in  auto-  and  cross- 
spectral  estimation  of  stationary  random  processes  by  means  of  overlapped  Fast 
Fourier  Transform  (FFT)  processing  are  derived  for  a  variety  of  constraints.  We  con¬ 
strain  the  time  duration  of  the  data  window,  or  constrain  the  bandwidth  in  various  ways, 
such  as  the  equivalent-noise  bandwidth,  the  half-power  bandwidth,  or  the  root-mean-square 
bandwidth.  In  the  three  bandwidth-constrained  cases,  the  window  duration  is  adiusted  to 
meed  the  constraint. 

We  find  that  the  Hanning  window  is  a  reasonable  compromise  for  achieving  minimum 
bias,  because  in  addition  to  being  the  optimum  for  one  bandwidth  constraint,  it  is  very 
close  to  the  optima  for  two  other  bandwidth  constraints.  The  relative  merits  of  the  spec¬ 
tral  characteristics  of  the  windows  are  also  discussed. 
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ABSTRACT 


The  time- limited  nonnegative  data  windows  that  minimize  the 
bias  in  auto-  and  cross-spectral  estimation  of  stationary  random 
processes  by  means  of  overlapped  Fast  Fourier  Transform  (F FT) 
processing  are  derived  for  a  variety  of  constraints.  When  the  time 
duration  L  of  the  data  window  is  constrained,  the  optimum  data 
window  is  (2/L)1//2  cos  (tt  t/L),  |t|  S  L/2;  when  the  equivalent- 
noise  bandwidth  is  constrained,  the  optimum  data  window  is 
(8/3L)1/2  cos^  (rrt/L),  which  is  the  Hanning  window;  when  the  half¬ 
power  bandwidth  is  constrained,  the  optimum  data  window  is 
L"1' 2  [1.682  +  4.261  cos  (4. 434  t/L)  -4.  337  cos  (3.552  t/L)], 
which  is  very  similar  to  the  Hanning  window;  and  when  the  root- 
mean-square  bandwidth  is  constrained,  the  optimum  data  window  is 
4/(5L)l/2  cos0  (it t/L).  In  the  three  bandwidth-constrained  cases, 
1’ie  window  duration  L  is  adjusted  to  meet  the  constraint. 

The  Hanning  window  is  a  reasonable  compromise  for  achieving 
minimum  bias,  because  in  addition  to  being  the  optimum  for  one 
bandwidth  constraint,  it  is  very  close  to  the  optima  for  two  other 
bandwidth  constraints.  The  relative  merits  of  the  spectral  charac¬ 
teristics  of  the  windows  are  also  discussed. 
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MINIMUM- BIAS  WINDOWS  FOR 
SPECTRAL  ESTIMATION  BY  MEANS  OF 
OVERLAPPED  FAST  FOURIER  TRANSFORM  PROCESSING 


INTRODUCTION 

The  selection  of  good  data  windows  in  spectral  estimation  of  stationary 
random  processes,  to  minimize  leakage,  is  an  important  consideration  and  has 
received  much  attention  [1-8].  In  [6],  a  thorough  investigation  of  fourgood  data 
windows  revealed  virtually  the  same  variance-reduction  capabilities  of  over¬ 
lapped  Fast  Fourier  Transform  (FFT)  processing  when  the  proper  overlap  was 
used  for  each  window.  The  ultimate  variance  reduction  of  this  direct  procedure 
was  also  demonstrated  to  be  identical  to  that  attained  by  the  older  (indirect) 
analysis  procedure  in  [4], 

In  this  report,  attention  is  focused  on  the  bias  in  the  estimation  of  power 
density  spectra  by  means  of  overlapped  FFT  processing.  Specifically,  the  bias 
is  minimized  by  the  choice  of  data  windows  that  are  restricted  to  be  time-limited 
and  nonnegative  and  are  subject  to  either  a  time-duration  constraint  or  a  band¬ 
width  constraint.  These  results  complement  and  extend  those  of  [8]  for  the  in¬ 
direct  approach  to  spectral  estimation. 


PROBLEM  DEFINITION 

The  overlapped  FFT  method  for  spectral  estimation  and  the  reasons  for  its 
use  are  documented  in  [6].  The  mean  of  the  spectral  estimate  is  given  by  [6, 
eq.  (5)]  for  auto-spectral  estimation,  and  by  [7,  eq.  (4A)]  for  cross-spectral 
estimation.  In  both  cases,  the  mean  takes  the  form* 

E  (G(f)}  =  /dv  G(f-v)  |  W(v) |2  ,  (1) 

A 

where  G(f)  is  the  estimate  of  the  true  (auto  or  cross)  spectrum  G(f),  and 

W(f)  =  fdt  exp(-i2irft)  w(t),  (2) 

♦Integrals  without  limits  are  over  the  range  of  nonzero  integrand. 


1 


TR  4513 


where  w(t)  is  the  data  window  multiplied  in  the  time  domain  by  the  available 
data.  It  is  assumed  that  the  data  window  is  time-limited  and  nonnegative: 


w(t) 


|  =  0  for  |t|  >L/2  1 
0  for  |t|  <L/2  r 


(3) 


where  L  is  the  window  duration.  The  restriction  to  nonnegative  data  windows 
guarantees  that  the  spectral  window  W(f)  peaks  at  the  origin.  It  should  be  noted 
that  (1)  is  true  with  no  restriction  on  the  available  record  length  T  and  with  no 
restriction  on  the  statistics  of  the  random  processes  involved,  except  that  the 
processes  must  be  stationary;  they  need  not  be  Gaussian  for  (1)  to  apply. 

The  desired  value  of  (1)  is  the  true  value  G(f);  therefore  the  bias  in  estima¬ 
tion  is  defined  as 


B(f)  =  E  { G(f) }  -  G(f).  (4) 

We  approximate  this  bias  by  expanding  G(f-i>)  in  (1)  according  to 

G(f-v)  S  G(f)  -vG'(f)  +  1/2  G"<f)v2  -  1/6  G'"(f)v3  +  1/24  G""(f)v4,  (5) 

where  the  prime  denotes  a  derivative.  Substitution  of  (5)  and  (1)  into  (4)  yields 
B(f)  S  1/2  G"(f) Jdv  v2 1  W(v)  |2  +  1/24  G""(f) fdv  v4|W(v)|  2,  (6) 

where  we  have  assumed  (without  loss  of  generality)  that 

/ dv  | W(v)  | 2  =/dt  w2(t)  =  1,  (7) 


and  that  |w(v)|2is  even  about  the  origin;  that  is,  w(t)  is  a  unit-energy  real 
waveform . 

We  express  (6)  as 

B(f)  S  1/2  G"(f)  DX  +  1/24  G""(f)  D2  =  Bi  (f)  +  B2  (f),  (8) 

where  window  constants 


D1  =7dvv2  | W(v)|2, 
D2=/dv/  |W(v)|2, 
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are  independent  of  the  true  spectrum  G(f).  To  minimize  the  bias,  we  must 
therefore  minimize  and/or  D2,  subject  to  (3)  and  (7).  To  this  aim,  it  is 
useful  to  express  (9)  in  terms  of  the  time  domain.  There  follows,  by  use  of  (2), 


and 


Dj  =  (2ir)-2  fdt  [w’(t)]2  (10) 

D2  =  (2ir)-4/dt  [w"(t)]2.  (11) 


Strictly,  the  approximation  (8)  to  the  bias  is  due  only  to  local  variations  in 
true  spectrum  G(f)  about  the  frequency  point  f  of  interest;  equation  (5)  is  not 
necessarily  a  good  approximation  for  larger  v.  Thus,  peaks  in  tbr:  true  spectrum 
that  are  distant  from  the  point  f  under  investigation  are  not  accounted  for  by  (5). 
To  minimize  the  effects  of  remote  soectral  peaks  on  bias,  we  must  also  require 
that  the  spectral  window  W(f)  decay  sufficiently  rapidly  for  large  |f|.  Thus 
the  results  of  the  following  optimizations  are  not  final,  but  must  be  investigated 
to  see  if  they  also  meet  the  requirement  of  sufficiently  rapid  decay  with  frequency. 

In  addition  to  constraints  (3)  and  (7),  we  shall  be  interested  in  constraining 
the  bandwidth  of  the  window;  this  is  in  keeping  with  the  philosophy  of  requesting 
a  specified  frequency  resolution  for  spectral  estimation,  and  letting  the  window 
duration  L  and  overlap  be  whatever  is  necessary  to  meet  this  requirement  [6], 

It  should  also  be  noted  that  constraints  on  bandwidth  tend  to  equalize  the 
variance-reduction  capabilities  of  the  windows.  This  may  be  seen  from  [6  , 
eq.  (22)],  where  the  equivalent  number  of  degrees  of  freedom  is  given  approxi¬ 
mately  by 


2TBst  for  T»L  ,  (12) 

where  Bst  is  the  statistical  bandwidth  [9,  p.  278]  of  the  window.  Thus,  if  all 
windows  were  constrained  to  have  the  same  statistical  bandwidth,  they  would  all 
have  the  same  variance-reduction  capabilities,  and  we  could  minimize  the  bias 
subject  to  this  constraint.  However,  this  constraint  is  not  mathematically  trac¬ 
table.  *  Therefore,  we  resort  toconstraints  on  other,  more  tractable,  bandwidth 
measures,  with  confidence  that  they  too  will  yield  comparable  variance  in  spec¬ 
tral  estimation  (see  [6,  table  1]). 

*We  have  not  been  able  to  express  the  time-domain  constraint  (3)  directly  in  the 
frequency  domain,  nor  have  we  been  able  to  express  the  requirement  that 
4»W(T)  be  a  legal  correlation  function  directly  in  the  time  domain;  see  [6, 
eqs.  (7)  and  (17)-(2l)]. 
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PROBLEM  SOLUTION 

Four  different  constrained  problems  will  be  addressed  in  this  section:  con¬ 
strained  window  duration  L  ;  constrained  equivalent-noise  bandwidth;  constrained 
half-power  bandwidthjand  constrained  root-mean-square (rms)bandwidth.  The  win¬ 
dow  duration  L  is  adjusted  to  meet  the  bandwidth  constraint  in  the  latter  three  cases. 

DURATION  CONSTRAINT 

Here  we  wish  to  minimize  Bj  in  (10),  subject  to  constraints  (3)  and  (7)  and 
a  fixed  value  of  window  duration  L.  In  order  that  (10)  be  finite,  w(t)  must  be 
continuous;  therefore  w(±L/2)  =  0  from  (3).  When  we  use  a  calculus-of- variations 
approach,  the  optimum  window  w0(t)  must  satisfy  the  differential  equation 

w<J(t)  +  Xw0(t)  =  0,  |t|  <  L/2,  (13) 

where  K  is  a  constant  (Lagrange  multiplier).  The  solution  of  (13)  that  satisfies 
the  boundary  conditions  and  (7),  and  has  minimum  D^,  is 

wD(t)  =  (-f;)1^2  cos(irt/L),  |t|  S  L/2.  (14) 

The  corresponding  value  of  (10)  is 


Several  windows  are  compared  in  table  1.  It  is  seen  that  the  Hanning  window 
has  33  percent  greater  bias,  as  measured  by  B^f),  than  the  optimum  window, 
urder  a  Juration  constraint. 


Table  1.  Window  Bias  Constants 


Data  Window 

4L2Di 

Optimum,  (14) 

Parabola 

Triangle 

Hanning 

1 

10/,r2  =  1.01 

12/tt2  =  1.22 
4/3  =  1.33 
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The  spectral  window  corresponding  to  the  optimum  data  window,  (14),  is 

Wo(f)^(2L)V2^il^l.  ,16) 

1-4L  i2 

o 

The  decay  for  large  frequencies  is  only  as  f  .  The  sidelobes  of  this  and  the 
following  windows  will  be  discussed  later. 

EQ  UIV  A  LENT-  NOISE-  BANDWIDTH  C  ONSTRAINT 


The  equivalent-noise  bandwidth  Be  of  spectral  window  W(f)  is  defined  as 

JfJwffilL.  !_  ,  (1?) 

|  W(0)  | 2  [7dtw(t)]2 

where  we  have  used  (7)  and  (2).  The  quantity  Bg  can  be  interpreted  physically 
as  the  bandwidth  of  an  ideal  rectangular  filter  that  would  pass  the  same  amount 
of  power  as  a  filter  W(f) ,  when  subjected  to  white  noise;  see  figure  1.  ThepeaK 
of  |W(f)|2  occurs  at  the  origin,  since  data  window  w(t)  is  nonnegative. 


Figure  1.  Equivalent-Noise- Bandwidth  Interpretation 


The  problem  here  is  to  minimize  in  (10),  subject  to  constraints  (3),  (7), 
and  (17).  TMs  problem  is  solved  in  appendix  A,  with  the  result  that  the  optimum 
data  window  is 


cos2(irt/L),  |t|  S  L/2  . 


(18) 
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This  is  the  familiar  Hanning  window.  The  window  duration  L  must  be  chosen  as 


L 


3 

2  Be 


9 


according  to  (17).  The  minimum  value  of  D-^  is 


(19) 


Di  = 


(20) 


The  specified  equivalent-noise  bandwidth  dictates  the  window  duration  L  and 
the  minimum  attainable  bias  constant  D^.  The  spectral  window  corresponding 
to  (18)  is  . 

I  2L  \  1/2 


w0(f)  -  - 


V  3 


sin(TrLf) 

irLf(l-L2f2) 


(21) 


where  L  must  be  determined  from  (19).  The  decay  for  large  frequencies  varies 
as  f  “3. 


HALF- POWER- BANDWIDTH  CONSTRAINT 


The  half-power  bandwidth  B^  of  spectral  window  W(f)  is  defined  as 


W(tBh/2) 

W(0) 


2 


(22) 


We  desire  to  minimize  Dj  in  (10),  subject  to  constraints  (3),  (7),  and  (22). 
Converting  (22)  into  the  time  domain  and  restricting  w(t)  to  be  even,  *  constraint 
(22)  takes  a  desirable  integral  form: 

fit  w(t)[cos(7rBht)  -  2-1/2]  =  0  .  (23) 

The  solution  to  this  minimization  problem  is  presented  in  appendix  B.  The  optimum 
data  window  is 


w0(t)  =  L~1^2  [1.682  +  4.261  cos(4.434t/L) 

-  4.337  cos(3.552t/L)],  |t|  S  L/2  .  (24) 


*  An  odd  component  in  w(t)  increases  the  rate  of  variation  and  therefore 
increases  Dj. 
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The  window  duration  L  must  be  chosen  as 

1.411 


according  to  (23).  The  minimum  value  of  Dj  is 

D]  =  0.1604  . 


(25) 


(26) 


The  spectral  window  corresponding  to  (24)  is  obtained  by  employing  (2);  this  will 
be  discussed  in  the  next  section.  It  is  shown  in  appendix  B  that  w^(iL/2)  =  0; 
therefore  the  decay  of  the  spectral  window  is  according  to  f"3. 

For  comparison,  if  the  time  duration  of  the  Hanning  window  is  adjusted  to 
realize  the  specified  half-power  bandwidth,  namely  L  =  1.441/B^  [6,  eq.  (33) 
and  table  l],  it  follows  that  Dj  =  0.1606  B^ .  Thus  the  Hanning  window  has 
virtually  the  same  bias  as  the  optimum  window  under  a  half-power-bandwidth 
constraint.  Further  comparisons  are  made  in  the  next  section. 

R  OOT-  M  EA  N-  SQUAR  E-  BA  NDWIDTH  C  ONSTRA INT 

The  rms  bandwidth  Br  of  spectral  window  W(f)  is  defined  as 


Inspection  of  (7)  and  (9)  immediately  reveals  that 

Dj  =  B^  .  (28) 

Thus  if  the  rms  bandwidth  is  constrained,  bias  constant  Dj  is  fixed.  In  this 
case,  it  is  reasonable  to  resort  to  minimization  of  the  second  bias  constant  D2 
in  (9)  or  (11).  Thus,  we  wish  to  minimize  (11),  subject  to  constraints  (3),  (7), 
and 

j/dt[w'(t)]2  =  (27rBr)2  .  (29) 

The  solution  to  this  problem  is  presented  in  appendix  C.  The  optimum  data  win¬ 
dow  is 

w0(t)  - - -  cos3(irt/L),  |t|  5  L/2  .  (30) 

(5L)1/2 
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The  window  duration  L  must  be  chosen  as 

3  1 


L  = 


2*75  Bt 


according  to  (29).  The  minimum  value  of  D2  is 

£,91 

The  so^ctral  window  corresponding  to  (30)  is 

W0(f)  =  21  (5L)1/2 _ ±°aJ?¥J _ 

15tt  (l-4L2f2)(l-4L2f2/9) 


(31) 


/QO\ 


(33) 


where  L  is  determined  from  (31).  The  decay  for  large  frequencies  is  according 
to  f"4. 


For  comparison,  let  the  time  duration  of  the  Hanning  window  be  adjusted  to 
realize  the  specified  rms  bandwidth  Bf.  Then  employing  (18)  in  (29),  we  find 
L  =  l/(*s/3  Br),  and  (11)  yields  D2  =  3  B4  .  Thus  the  Hanning  window  has 
8  percent  more  bias  than  the  optimum  window  under  an  rms  bandwidth  constraint, 
as  measured  by  bias  constant  D2. 


COMPARISON  OF  WINDOW  CHARACTERISTICS 

In  figure  2,  one- half  of  the  symmetric  optimum  data  windows  for  the  three 
bandwidth-constrained  cases  are  drawn  for  a  common  time  duration  of  L  =  1. 
The  equivalent-noise-bandwidth  data  window  (Harming)  and  the  half-power-band¬ 
width  data  window  are  virtually  identical  and  are  continuous  in  value  and  deriv¬ 
ative  at  0.5.  The  rms-bandwidth  data  window  is  more  peaked,  and  goes  to  zero 
in  value,  in  derivative,  and  in  second  derivative  at  0.  5.  Thus  the  last  window 
would  require  greater  overlap  than  the  first  two,  in  order  to  realize  the  same 
variance  reduction;  see  [6]. 

In  order  to  deduce  the  required  overlap  for  the  rms  bandwidth  data  window, 
the  quadratic  and  cubic  data  windows  [6,  pp.  10-18]  are  superposed  in  figure  3. 
Over  most  of  the  range,  the  quadratic  and  rms-bandwidth  windows  are  very  close. 
Near  the  end  of  the  range,  however,  the  taper  of  the  rms-bandwidth  window  ap¬ 
proaches  that  of  the  cubic;  in  fact,  both  are  continuous  in  second  derivative  at 
0.5.  Thus,  it  is  anticipated  from  earlier  results  [6,  table  4  ]  that  slightly  over 
65  percent  overlap  would  be  required  for  the  rms  bandwidth  data  window  to  realize 
99  percent  of  its  maximum  equivalent  degrees  of  freedom  [  4,  p.  22  ]. 
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2. Or 


•  MS 
EN 


l>) 


•MS:  ROOT -MEAN-SQUARE-BANDWIDTH  DATA  WINDOW 

EN  EQUIVALENT-NOISE-BANDWIDTH  DATA  WINDOW 

HP  HALF-POWER-BANDWIDTH  DATA  WINDOW 


\ 


V 


•MS  EN 


HP 


0.5 


Figure  2.  Three  Bandwidth-Constrained  Data  Windows; 
L  =  1,  Unit  Energy 


o.o  ,  0.3 


Figure  3.  Comparison  of  RMS-Bandwidth  Data  Window  with 
Quadratic  and  Cubic  Data  Windows;  L  =  1,  Unit  Energy 
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The  spectral  characteristics  of  the  four  windows  derived  in  this  report  are 
presented  in  figures  4  through  7.  The  abscissas  on  every  plot  are  in  units  of  the 
half-power  bandwidth;  thus  all  the  curves  go  through  half-power  (-3.01  dB)  at 
f/B^  =  0.5.  The  duration-limited  window,  (14),  is  plotted  in  figure  4  and  ex¬ 
hibits  relatively  slow  decay  with  frequency,  since  the  data  window  is  discontin¬ 
uous  in  derivative  at  its  edge.  The  equivalent-noise-bandwidth  (Hanning)  and 
half-power-bandwidth  spectral  windows,  plotted  in  figures  5  and6,  are  virtually 
identical  and  have  good  decay  with  frequency,  since  the  data  windows  are  con¬ 
tinuous  in  derivative  at  their  edges.  The  spectral  window  for  the  rms-bandwidth 
data  window  is  plotted  in  figure  7  and  exhibits  very  rapid  decay  with  frequency. 
However,  as  noted  above,  by  virtue  of  requiring  greater  overlap  for  maximum 
variance  reduction,  this  window  will  require  somewhat  greater-size  FFTs  than 
do  the  other  windows. 


CONCLUSIONS 

The  Hanning  window  is  optimum  under  an  equivalent-noise-bandwidth  con¬ 
straint,  as  far  as  minimization  of  bias  constant  Dj  is  concerned.  Furthermore, 
it  is  near  the  optima  for  two  other  bandwidth  constraints.  Its  spectral  decay  is 
also  sufficient  for  most  cases  that  the  bias  is  relatively  unaffected  by  distant 
spectral  peaks.  And  with  50  percent  overlap,  it  realizes  92  percent  of  the  max¬ 
imum  number  of  equivalent  degrees  of  freedom  [6,  table  6].  Thus,  the  Hanning 
window  is  a  reasonable  compromise  to  utilize  in  spectral  estimation  of  random 
stationary  data. 
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Figure  5.  Spectral  Window  for  Half- Power- Bandwidth  Data  Window 
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Figure  7.  Spectral  Window  for  RMS- Bandwidth  Data  Window 
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Appendix  A 

DERIVATION  OF  THE  OPTIMUM  DATA  WINDOW  FOR  AN 
EQUIVALENT-NOISE-BANDWIDTH  CONSTRAINT 

We  wish  to  minimize  (10)  (from  the  main  text),  subject  to  constraints  (3), 
(7),  and  (17).  Equations  (7)  and  (17)  are  integral  constraints  and  are  in  a  con¬ 
venient  form  for  a  calculus- of- variations  approach.  The  way  we  handle  (3)  is  to 
first  ignore  the  nonnegative  limitation;  then,  out  of  the  class  of  allowable  sol¬ 
utions,  we  restrict  attention  only  to  the  nonnegative  solutions  and  pick  the  best. 

In  order  that  (10)  be  finite,  w(t)  must  be  continuous.  Using  (3),  we  see  that 
this  means  that 


w(±L/2)  =  0.  (A-l) 

A  calculus-of-variations  approach  tells  us  to  minimize  the  quantity 

Q  =  JjtjV (fiT  -  \  [#  «/ft)  <-  2\ Wit),  (A-2) 

where  \  and  are  Lagrange  multipliers;  the  resulting  differential  equation 
for  the  optimum  window  is 

*>)  +  =  V,  W<  L/2.  <A-3> 


We  employed  (A-l)  on  the  allowed  variations  in  deriving  (A-3). 
The  general  solution  of  (A-3)  is 


A  cos(at)+  B  Jin  (#■£)  4-  C 


OR 


\t)  A  coal*  (of)  4-  B  5inh(ert)  +  C 

OR 

A  +  Bt  J 


>  Krl<  l/2  . 


(A- 4) 


where  a  is  real  and  positive.  The  third  alternative  in  (A-4)  yields  the  trivial 
solution  when  (A-l)  is  imposed. 
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The  second  alternative  reduces  to 

=  A[cosl»(<rt)  -  cos  >•(*0)],  W  s  L/2' 


(A-5) 


Since  A  can  be  chosen  to  satisfy  the  energy  constraint  (7),  and  L  can  be  chosen 
to  satisfy  the  bandwidth  constraint  (17),  there  is  left  only  the  variable  d  in 
(A-5)  to  vary;  ft)  is  certainly  nonnegative  by  the  choice  of  proper  polarity 
for  A.  In  order  to  find  the  best  value  of  ft  for  this  second  alternative  of  (A-4), 
we  compute  P,  in  (10)  versus  A  and  pick  the  minimum. 

To  accomplish  this  goal,  we  define 

JjpU)  =  !=  ^  |«1  S  I,  (A-6) 


where 


and 


u)  =■  cosh  W“)  ~  coshM }  M ~  1 


*  -  aL/2. 

The  bandwidth  constraint  (17)  then  becomes 


where 


K#  S  f  Ji  (a). 


-  l 

The  energy  constraint  (7)  becomes 

a  u 


where 


A  2^*  i, 

k  ?  r'dui‘(«). 


The  window  constant  (10)  becomes 

T)  «  A.  u' 

^  2tt‘  L  > 


(A- 7) 

(A-8) 

(A- 9) 

(A- 10) 
(A— 11) 

(A-12) 

(A-13) 
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where 

5  du  [A7h)]  . 


(A- 14) 


If  we  eliminate  A  and  L  by  using  (A-9)  and  (A-ll),  then  (A-13I  becomes 


T>, 


KJ  •* 


(A- 15) 


The  quantities  necessary  in  (A-15)  follow  from  (A-7),  (A-10),  (A-12),  and  (A-14): 

K.=  2[c6sV,W-^I]) 

is  -  .  il 

n  ~  j  ^  (A-i6) 

K..  2 


When  (A-16)  is  substituted  in  (A-15),  T>,  is  found  to  increase  monotonically  with 
increasing  4  ;  the  limit  as  o(-*  H  is 


US'  V2  _ 
12tt*  Be 


,nn 


(A-17) 


The  first  alternative  in  (A-4),  when  subjected  to  boundary  condition  (A-l), 
breaks  into  two  subcases.  In  the  first,  if  Sin(*)  *  0,  then  "B  is  arbitra-y. 
Then  K>  I  .  The  function  with  smallest  D,  corresponds  to  K*t, 

and  yields 

W,  It)  -  A[c0S (>wt/L)  +  l]  +  8  9'n (2irt/t.),  1*1  *  L/ 2  *  (A"18) 

However,  we  must  have  B  =  0  ;  otherwise  W0(t)  would  go  negative  somewhere. 
Imposition  of  the  energy  constraint  (7)  and  bandwidth  constraint  (17)  yields 

W,(t)  COS*(Tt/L),  HI  i  0>  (A_19) 

where 


L 


(A- 20) 
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The  corresponding  value  of  Dj  is 

which  is  smaller  than  (A-17). 


(A-21) 


The  second  subcase  occurs  if  Sinfr)  *  0  •  Then  B  in  (A- 4)  must  be 
zero,  and  we  get 


W0  (t)  *  “  C0S(*3,  W  5  U2> 


(A- 22) 


where  •<  k  IT  .  The  comments  and  method  immediately  below  (A-5)  apply 
equally  well  here.  Therefore  we  define 


-  cos^u)  -  cos("^;  |u|  i  I,  <  ^  Ktt, 


and  find 


_  f  s'in  itf) 

‘‘I  * 


K*  -  2 


2  2*C  J  * 


(A-23) 


(A- 24) 


When  (A-24)  is  substituted  in  (A-15),  J>,  is  found  to  decrease  monotonically  with 
increasing  •*;  ,  at  least  for  oi  up  to  tt.  However,  when  oOTT,  -4(u)  becomes 
negative  somewhere;  this  may  be  seen  by  noting  that  =  -°<  sin  (p<)  is 

positive  if  <>■"“.  Thus  the  limiting  member  of  this  subclass,  which  is  (A-19),  is 
the  outimum  window. 
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Appendix  B 

DERIVATION  OF  THE  OPTIMUM  DATA  WINDOW  FOR  A 
HALF- POWER-BANDWIDTH  CONSTRAINT 


Here  we  will  minimize  (10),  subject  to  constraints  (3),  (7),  and  (23).  In 
order  that  (10)  be  finite,  (A-l)  must  again  be  true.  A  calculus- of- variations 
approach  tells  us  to  minimize 

Q  =  (b-1) 

where  ^ 

c  (t)  s  cos^t)  -(ryW5'TrB^)®,=2  ,  (B-2) 

and  Xj  and  Xa  are  Lagrange  multipliers;  the  resulting  differential  equation  for 
the  optimum  window  is 

W,"  ttj  +  X,  W.i-fc)  *  X2  cf*-\  Hr)  <  (B-3) 


If  X,  is  negative,  the  form  for  wcfc)  includes  sinh  and  cosh  terms,  which 
lead  to  a  progressively  larger  value  of  D,  as  X,  becomes  more  negative, 
similar  to  the  result  of  appendix  A.  If  X,  is  zero,  the  only  solution  to  (B-3) 
and  (A-l)  is  the  trivial  solution.  If  X,  is  positive,  the  general  solution  to  (B-3) 
is 


Klt)  = 


OR 

A  cos  (u>t)  +  Bain  (<o  *-) ;  X  =  w 


(B-4) 


We  discard  the  odd  solutions  for  the  reason  given  in  the  footnote  to  (23). 

If  we  attempt  to  use  the  second  alternative  in  (B-4),  we  can  eliminate  A 
by  means  of  energy  constraint  (7).  However,  the  bandwidth  constraint  (23)  can 
not  be  satisfied  for  any  value  of  oo  L  f  0  •  Therefore,  we  must  discard  the 

second  alternative. 

To  handle  the  fir  alternative  in  (B-4)  conveniently  ,,  we  define  a  function 

3o  (m)  s  W#  (^-u)  =A  cos(<m)  “  y  & ,  1^! < 
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where 


•BfcL. 


(B-6) 


Imposition  of  (A-l)  forces  (B-5)  into  the  form 


where  C  is  a  constant.  The  energy  constraint  (7)  yields 

c  --  (iVf  >  (B-8) 

where  we  have  used  (B-5)  and  (B-7),  and  defined 

K,  =  \  du  A  (u) .  (B-9) 

Satisfaction  of  bandwidth  constraint  (23)  demands  that 

0  -  j  du  [cos^m)  -  (B_10) 

where  we  employed  (B-6)  and  (B-7).  Substitution  of  the  detailed  form  for  -Afv), 
(B-7),  into  (B-10)  yields  the  relation 

=  j^cosk) +&(»-£)]  _5(*V*)  +’  5(^-  (B"U) 


(B-8) 


(B-9) 


where 


5(y)  * 


Sln(X) 


(B-12) 


For  a  given  value  of  (B-ll)  must  be  solved  for  the  smallest  value  of  *  ;  £ 
is  a  known  specified  constant.  In  order  to  find  the  best  value  of  ,  we  com¬ 
pute  D,  versus  and  pick  the  minimum,  always  being  careful  that  Jtf^remain 
nonnegative  for  all  lull  l .  The  quantity 

t>» 


1  JKjl. 

VL*  K# 


Ki  p 1 

4S\  Bh 


(B-13) 
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using  (B-5),  (B-7),  (B-8),  (B-6),  and  defining 

K,  =  JduCVMl . 


(B-14) 


The  quantities  K.  and  K,  are  available  upon  substitution  of  (B-7)  in  (B-9)  and 
(B-14).  There  follows,  upon  use  of  (23), 

k,  =  c‘[i4  sivi)]  ♦ 

+  c,  c,  [2  a-  (2  -  f )  5  kt) + S  K+  -0  +  SH"  *3 ,  ,B"15) 

81  -  +  v1  c;ji-  s(*>] 

4-  2^C,  C^sl^-  «)  -  §(•<■),+  -)!  J ,  (B-16) 

where 

Cl  *  fosi^pT'  ,  C2  =  -  [(^cosk)  4-  Gr ()-£)]  .  (B— 17) 

The  numerical  approach  may  now  be  summarized  as  follows.  A  value  of 
0^  is  picked,  and  (B-ll)  is  solved  for  ot  .  Then  (B-17)  is  computed  and  sub¬ 
stituted  in  (B-15)  and  (B-16),  thereby  enabling  evaluation  of  D,  in  (B-13)  for 
that  choice  of  0  .  (Up  to  this  point,  6r  could  be  any  desired  constant;  we  now 
restrict  6r  ■  ).  When  this  approach  is  tried,  it  is  found  that  P, 

increases  monotonically  with  increasing  <^  .  On  the  other  hand,  when  is 
made  too  small,  becomes  negative  near  ia*|.  The  optimum  value  of  ^ 
is  realized  when  >■(■) -0.  From  (B-7)  and  (B-17),  this  requirement  is 


C,  5m 4-  C,  ^  Sink)  -  0. 


(B-18) 


The  simultaneous  solution  of  (B-ll)  and  (B-18),  with  smallest  *  ,  is  then  given 
by 

<p. ion  1116  ,  <  =  2.217  osir, 

The  optimum  value  of  P,  then  follows  from  (B-13)  as 

P,  ^  .\6CH$J0  K  , 

and  the  segment  length  follows  from  (B-6)  as 


1  -  Al.  J _  = 

L  tt  Bh 


Mil  423? 


Finally,  the  optimum  window  follows  from  (B-5)  as 

W#W-  67j|  4-1.261  0611  C^(4.13l!l1ft/0 

-  1, 337  3*91  C05(3. 552  26 13  t/0] ,  lii  *  L/2. 


(B-19) 


(B-20) 


(B-21) 


(B-22) 


21/22 
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Appendix  C 

DERIVATION  OF  THE  OPTIMUM  DATA  WINDOW  FOR  AN  RMS- 
13  ANDWIDTH  CONSTRAINT 


The  problem  here  is  to  minimize  (11),  subject  to  constraints  (3),  (7),  and 
(29).  In  order  that  (11)  be  finite,  it'lb)  must  be  continuous.  Using  (3),  this  means 
that 

w(±L/l)  =■  0,  w'(iL/l)=0,  (C-l) 

A  cal  cuius- of- variations  approach  tells  us  to  minimize  the  quantity 

Q  X^d-t^w'K:)]  +- (c“2) 

where  X  and  ^  are  Lagrange  multipliers;  the  resulting  differential  equation  for 
the  optimum  window  is 

•Tip  -  XK"(t)  ♦rw.lti  =0,  Itl  <  L/2.  (C-3) 

We  employed  (C-l)  on  the  allowed  variations  in  deriving  (C-3). 


To  solve  (C-3),  we  assume  a  form  e*p(st)  for  Substitution  in  (C-3) 

requires  that  s  be  chosen  to  satisfy 

S4-  XsV  l>  =  0  ;  i-  =  t  [X  ±  (X*-  fyf] .  (C-4) 


*>  "  z. 


(C-4) 


At  this  point,  several  alternatives  are  possible.  The  first  case  we  pursue 
is  a  negative  discriminant  in  (C-4).  Then  the  four  values  for  3  in  (C-4)  can  be 
expressed  as  *  * 

S  =  2,2  /-*,  >  (C-5) 

where  i  is  a  complex  constant  with  nonzero  imaginary  part.  For  distinct  roots 
(i.  e. ,  *  not  purely  real),  the  optimum  window  is 

w,  it)  -  f\  e*p(*t)  +-B  t*p(z-* i)  4-  C  e  xp ('*£■) + D  **p("2* t).  (C-6) 


In  order  that  (C-l)  be  satisfied  with  a  nontrivial  solution,  it  is  necessary  that 
the  determinant 

E  c  '/f  '/e* 

'/e  Ye?  E  £* 

z£  -2/E  -£/«*  (c_. 

f/£ 
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5=  < 


equal  zero,  where  E  s  exp(fcl-/2)  .  This  requires  that 

-  °-  <c'8) 

But  none  of  the  factors  in  (08)  can  be  zero  without  double  roots  resulting  for  s, 
in  which  case  the  form  (06)  is  not  appropriate.  Therefore,  the  negative- 
discriminant  case  is  self-destroying. 

The  second  case  corresponds  to  a  zero  discriminant  in  (04).  Then  we  have 
three  subcases: 

a,  a, 

ia,  iV-ia,  -iai  ,  «  real  and  positive.  (C-9) 

OK 

0,  0,0 ,  0  w 

For  subcase  1,  the  form  for  W#|±)  is 

H  It j  =  f\  cosk(«t)  4-  B  4- C  t  Cosli  (at)  4- Dt  sinK(at).  (c-10) 

Imposition  of  (C-l)  requires  that 

5inK{aL)=  ±  aL  (Oil) 

for  a  nontrivial  solution.  But  (Oil)  has  no  solutions  for  positive  real  q  . 

For  subcase  2  in  (09),  the  form  for  *(,10  is 

Hit)  =  A  to  s  (at)  4-B  sin  ( at)  4-  Ct  <®s(«t)  +Dt  -^(at).  (C-i2) 

Imposition  of  (C-l)  requires  that 

sin(al)B±aL  (C-13) 

for  a  nontrivial  solution.  Again  this  is  disallowed. 

For  subcase  3  in  (C-9),  the  form  is 

Hit)  =  A  4-  Bt  4  CtX4-  J>t3.  (C-14) 

Imposition  of  (C-l)  yields  only  the  trivial  solution. 

The  third,  and  last,  case  we  must  now  consider  is  a  positive  discriminant  in 
(C-4).  Then  we  have  three  subcases: 
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r  * «/  *  l~) 

3  ■  <  !fc  *  i  b  , 

I 


^  and  k  real , 


(C— 15) 


For  subcase  1,  we  have 


W,ffc)  =  A  (OsK(qt)  F  B  §ml»(ot)  +  C  COsMt't)  +  D  (C-16) 


Imposition  of  (C-l)  requires  that 


where 


. i  ,  ,„x  +a*h(«)  _ 

*  f  wip)  °r  * —  f  * 

•(  5  AL/2,  f  5  bL/z. 


(C-17) 


(C-18) 


The  only  solutions  of  (C-17)  are  ^  ;  that  is,  A  -  b  .  However,  these  have 

been  considered  already  in  (C-10)  and  found  inadequate. 

For  subcase  2  in  (C-15),we  have 

W#  X  /{  +  B  5 inb (art)  F  C  +•  X>  s'w(i>t).  (C-19) 


Satisfaction  of  (C-l)  demands  that 


fonhU)  •fa>n(^) 

*  =  P 


(C-20) 


The  second  alternative  in  (C-20)  leads  to  odd  solutions  only,  in  (C-19),  and  they 
must  be  discarded  because  of  their  higher  variation  rate.  The  first  alternative 
in  (C-20)  leads  to 


w  u.)  -  A  _ jt i  <  ^ 

WflH"  4  <**{?>  J 


(C-21) 


which  is  a  legal  nontrivial  solution.  The  values  of  o<  and  ^  are  related  as  shown 
in  figure  C-l.  To  handle  this  alternative  conveniently,  we  define  a  function 

. /LA  J  “AM  uslpj  1  .  .  ,  _ 


3.'-) 


=  =  c  - 


cosM*<) 


j  5  Cjklu),  I-  (<>22> 
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Figure  C-l.  Relationship  of  a  and  (J  in  (C-20) 


Then 


£(«)  =  T a»  =(tJ 


The  energy  constraint  (7)  requires  that 

./Uf 


where 


C  *&)  ■ 

Kj  s  j*  Ju  [-Ay)C»o] 


(C-23) 

(C-24) 

(C-25) 


The  bandwidth  constraint  (29)  requires  that 


(C-26) 


where  we  have  employed  (C-23),  (C-24),  and  (C-25).  Then  we  use  (C-23)  through 
(C-26)  to  determine  the  bias  constant  T>x  as 


K,K, 


(C-27) 


For  the  current  example  in  (C-22),  we  evaluate 

U  -  *(.*>»  i  .  ilM., 

0  COS*^) 
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(C-28) 


(029) 


For  an  arbitrary  ,  we  solve  the  first  alternative  in  (020)  for  p  , 

and  then  compute  Da  ,  by  means  of  (027)  and  (028).  It  is  found  that  Da 
increases  monotonically  with  orf  .  The  minimum  is  realized  when  e(  *  o  . 
namely 


(030) 


The  third  subcase  in  (015)  yields  the  form 


W0  It}  =  A  Gos(o-fc)  +  B  5l‘n  Ivt)  4-  C  (pi +■!>  (^)  •  (031) 


The  boundary  conditions  (C-l)  demand  that 


-fow  (*)  *■  ^ 


"  =  f 


(032) 


However,  the  second  alternative  in  (032)  yields  odd  solutions  and  is  discarded. 
The  first  alternative  yields  _  . 

<0-3, 

The  values  of  •(  and  ^  are  related  as  shown  in  figure  02.  As  above. 


“J;  W  -2:  >  (033) 


Figure  02.  Relationship  of  a  and  0  in  (032) 
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we  define  a  function 


(C-34) 


We  can  now  employ  (C-25)  througli  (C-27)  immediately.  We  evaluate 

K  =  _2  skyhsb-f)  i  f  3ty>) 

*  cos’ll  cos K)  vs((k)  Gosfy) 


For  an  arbitrary  X  (  ^  tt/ z)  ,  we  solve  the  first  alternative  in  (C-32)  for 
A  ,  and  then  compute  J>2  using  (C-27)  and  (C-35).  It  is  found  that  J>2 
decreases  monotonically  with  increasing  ^  ,  at  least  for  up  to  ir/2  ,  the 
limit  being  25/9  at  ^  =  rr/z  -  .  However,  when  o(  >  tt/2  ,  JK«)  goes 

negative  somewhere  and  is  unacceptable.  For  o(  *  rr/2  ,  (C-32)  is  not  an 
adequate  form;  we  note  instead  that  3  =  3ir/2  from  (C-31),  and  then 


iW  =  A  c.o s (5 u)  -  3 cos(£r  u) . 

(C-36) 

The  boundary  conditions  (C-l)  force  A  -  -3B  ,  giving 

Jki 'U  ~  C053(¥w)- 

(C-37) 

We  then  find 

K0  « 

y  k  -  -kd  .  k,  =  15X. , 

~J~>  ~  32  >  r*a-  Iz? 

(C-38) 

yielding 

L  =  5JF  b7  ’ 
b;, 

(C-39) 

which  is  smaller  than  (C-30).  Thus  (C-37)  is  the  optimum  window.  (Notice  that 

-A’o)  =  o.) 
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